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We study the quantum phase transition between Abelian and non-Abelian phases in an extended 
Kitaev spin model on the honeycomb lattice, where the periodic boundary condition is applied by 
placing the lattice on a torus. Our analytical results show that this spin model exhibits a continuous 
quantum phase transition. Also, we reveal the relationship between bipartite entanglement and the 
ground-state energy. Our approach directly shows that both the entanglement and the ground-state 
energy can be used to characterize the topological quantum phase transition in the extended Kitaev 
spin model. 



I. INTRODUCTION 



studied for the Hamiltonian: 



H 



J \ * X X 



n,m u n+l,m 



Quantum phase transitions, which occur when a driv- 
ing parameter in the Hamiltonian of the system changes 
across a criticalpoint, play a central role in condensed 
matter physics 0, Q • While most quantum phase transi- 
tions can be characterized by symmetry breaking, there is 
also an exception that can only be witnessed by topologi- 
cal order (see, e.g., [3, 4]). Signatures of topological order 
in many-body quantum systems can characterize a topo- 
logical quantum phase transition and include, e.g., the 
existence of excitations obeying fractional statistics (see, 
e.g., @), ground-state degeneracy related to the topology 
of the system (instead of the symmetry Wsee, e.g., |3],@|), 
and topological entanglement entropy [3, Q . In particu- 
lar, the spectral Chern number Q serves as a topological 
number for characterizing a two-dimensional (2D) system 
of noninteracting (or weakly interacting) fermions with 
an energy gap. Without closing the gap, energy spectra 
with different Chern numbers cannot be deformed into 
each other jig ]. This is because a topological quantum 
phase transition occurs when changing the Chern num- 
ber. 



Recently, it was shown |ll|, Il2| that the topological 
quantum phase transition in the Kitaev spin model can 
be characterized by nonlocal-string order parameters. In 
an appropriate dual representation, this order parameter 
can become local and the basic concept of Landau theory 
of continuous phase transition is also applicable [ill ]. 
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where a* m , m and m are the Pauli matrices at the 
site (n, m), with column index n = 0, 1, 2, 3, • • • , N — 1 
and row index m = 0, 1, 2, 3, • • • , M — 1. A nice Jordan- 
Wigner transformation was introduced [HI, [HI, Ell to 
solve this model and the redundant gauge degrees of free- 
dom were removed. 

The phase diagram of the Kitaev model (p} consists of 
two phases: A band insulator phase and a topologically 
non- universal gapless phase 0]. The insulator phase, as 
Kitaev has shown by using perturbation theory [jj Ei| , 
is equivalent to a toric code model [l6(. While Abelian 
anyons can be defined in the insulator phase, the vortices 
in the gapless phase do not have a well-defined statis- 
tics. Applying an external magnetic field as a pertur- 
bation, which breaks the time-reversal symmetry in Eq. 
([T]), a gap opens in the gapless phase and the vortices 
then obey a well-defined non-Abelian anyonic statistics 
Q. The third-order perturbation corresponds to exactly 
soluble models 0, Ecfl whose spectrum has recently been 
extensively studied [l7| . 

In this paper, we study the following Hamiltonian [ToL 
MM: 



In the Kitaev model, a ^-spin is placed at each site of a 
honeycomb lattice [see Fig. Ola)] and the interactions be- 
tween nearest-neighbor spins are highly anisotropic with 
three types of bonds J x ,J y , and J z . To simplify the 
site-labelling of the honeycomb lattice, one can deform it 
to a topologically equivalent brick-wall lattice shown in 
Fig. Qlb) . In [ll|, EH, EH , the topological quantum phase 
transition of the Kitaev model on a brick-wall lattice was 
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Hereafter, we call the model in Eq. ([2]) an extended Ki- 
taev model. We solve this model on a torus and mainly 
focus on the quantum phase transition between the phase 
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FIG. 1: (Color online) (a) Honeycomb lattice constructed by 
two triangular lattices that are connected to each other by 
three types of bonds J x ,J y and J z . (b) The brick- wall lat- 
tice, which is deformed from the honeycomb lattice in (a). 
This deformed lattice can be used to label the sites of the 
honeycomb lattice by column and row indices. 
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with Abelian anyons and the phase with non-Abelian 
anyons. We first apply the Jordan- Wigner transforma- 
tion to the spin operators and then introduce Majorana 
fermions to get the ground state of Eq. ((2|) in the vortex- 
free sector. We show that the third directional deriva- 
tive of the ground-state energy is discontinuous at each 
point on the critical line separating the Abelian and non- 
Abelian phases, while its first and second directional 
derivatives are continuous at this point. This implies 
that the topological quantum phase transition is contin- 
uous in this extended Kitaev model. Moreover, at this 
critical point, we also study the nonanalyticity of the 
entanglement (i.e., the von Neumann entropy) between 
two nearest-neighbor spins and the rest of the spins in the 
system. We find that the second directional derivative of 
the von Neumann entropy is closely related to the third 
directional derivative of the ground-state energy and it 
is also discontinuous at the critical point. Our approach 
directly reveals that both the entanglement measure and 
the ground-state energy can be used to characterize the 
topological quantum phase transition in the extended Ki- 
taev model. 



II. TOPOLOGICAL QUANTUM PHASE 
TRANSITION 
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FIG. 2: (Color online) Phase diagram of the extended Ki- 
taev spin model, where J x , J y and J z > 0. The gray region 
corresponds to the non-Abelian phase and the three trian- 
gular (light gray) regions correspond to the Abelian phase. 
The thick solid, dashed and dotted lines are A x — 1 + A y , 
A y = 1 + Ax, and A y = 1 — A^, where A x = Jx/Jz and 
A y = Jy/Jz. These lines consist of the boundary of the gray 
region, which are the critical lines separating the Abelian and 
non-Abelian phases. The thin dotted line intersects the thick 
solid and dotted lines at the points (A x , A y ) = (1.5, 0.5) and 
(0.5, 0.5). The direction I has an inclination angle <p with 
respect to the horizontal axis and it indicates the direction 
along which the driving parameters A^ and A y vary. 
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for n + m equal to an even integer. When the phase 
(arising from the Jordan- Wigner transformation) related 
to each bond between the (N — l)th column and the 
zcroth column is chosen to be 2irl (I is an integer ), the 
Hamiltonian is reduced to 



Let us define the Jordan- Wigner transformation [l| 
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where s = 1 if the integer n + to is odd and s = 2 if the 
integer n + m is even. Also, we introduce the following 
definitions for Majorana fermions: 
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In Eq. ©, the ^NM operators idn)md^ m+1 , where n + 
to is an even integer, commute with each other. The 
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ground state is in the vortex-free sector [Tj, [2(| with 

^L+i^n%^i+2,m+i^n+2,m = — 1j which corresponds to 
the case with the eigenvalue of each plaquette operator [9( 
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equal to 1. Thus, we can set the \NM operators 

id^md^m+i all equal to 1 in Eq. ([6|), in order to ob- 
tain the ground-state energy. For this quadratic Hamil- 

tonian, the Fourier transformation of H via Cnln — 

2 'Ei 
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where a x , a y and a z are Pauli matrices, <E> 
with c°j> = (c£A , and 
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Let us define 
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where 



■ fk x I fay \ ■ ( k x k y \ . - ky 

t(k) = iJ x e\ 2 ^> + iJ y e \ 2 ^> - iJ z e V5 , 



and 



e(k) = |h(k) 



= y|a(k)| 2 + 4J 2 sm 2 k x 
It is straightforward to verify that 
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i.e., i?£ and i?k are fermionic operators, and the Hamil- 
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tonian ([8]) can be written as 

H = Y J [e(k) - 2e(k)i?ti3 k 



For Hamiltonian p4|) . the ground-state energy is 
— J2k e (k) and the ground-state \g) obeys BjB k \g) = . 
The energy spectrum e(k) is gapless [10| only when 
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corresponds to the thick solid, dashed, and dotted lines 
in Fig. [2J respectively. 

When J > 0, the spectral Chcrn number is 1 if 

Jr ^ Ji 



J z , Jy < Jz + J 2 and J 2 < J x 
if J x > J y + J z or J y > Jz + Jz or J z > J x + J y (see 



J y , and 



[t|[Tc|). These two cases correspond to the non-Abelian 
and Abelian phases in the Kitaev model and both of 
them are gapped topological phases. The phase diagram 
is shown in Fig. [2] where the gray area corresponds to 
the non-Abelian phase and the critical lines (denoted as 
thick solid, dashed and dotted lines) separate the Abelian 
and non-Abelian phases. This indicates that the system 
can experience quantum phase transitions across these 
three thick lines. Here we rescale the inter-spin coupling 
strengths by introducing A^ = J x /J z ,A y = J y /J z , and 
A = J/Jz, so as to conveniently characterize the quan- 
tum phase transition. 

To demonstrate the quantum phase transition, one 
may reveal the nonanalyticity of the ground-state energy. 
The ground-state energy per site is 
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where BZ denotes the first Brillouin zone. Its directional 
derivatives with respect to the driving parameter along 
any given direction I (see Fig. [2]) are 
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If the nth directional derivative d n E/dl n (n=l, 2, • ■ ■) 
is nonanalytical at the critical point (A x ,A y ), and the 
directional derivatives d m E/dl m with < m < n are 
analytical there, a topological quantum phase transition 
occurs at this critical point. 

It can be proved that (see Appendix A) 
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FIG. 3: The ground-state energy per site E (in units of J 2 ) 
and its first, second, and third derivatives with respect to A^, 
where A y = 0.5, A = 0.1 and <p — (which corresponds to 
the horizontal thin dotted line in Fig. [2}. It is clear that 
E and 8E/dA x are continuous functions, but d z E/dA x is 
discontinuous at the transition points A^ = 0.5 and 1.5. 



denotes (A x — A y ) — > 1 with (A^ - A y ) < 1. In (Qjj]) 
Qi = i(4A 2 + A2+A, + l), 
Q 2 = \yj (8A 2 + 2A 2 , + 2A y - if + 3(1 + 2A y ) 
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Equations (fT"7]) and (fl"8|) reveal that a continuous topo- 
logical quantum phase transition occurs across the criti- 
cal line A x = 1 + A y (denoted by the thick solid line in 
Fig. [2]). Similarly, it can be shown that such a continuous 
topological quantum phase transition also occurs across 
the critical lines A y = 1 + A^ and A x = 1 — A y (de- 
noted, respectively, by the thick dashed and dotted lines 
in Fig. [2]). As a numerical test, we choose A = 0.1, ip = 
and Ay = 0.5 to show this quantum phase transition 
in Fig. [3] where the range of A x is chosen by the thin 
dotted line in Fig. It can be seen in Fig. [3] that the 
ground-state energy and its first and second directional 
derivatives are continuous for each A x , while its third di- 
rectional derivative is nonanalytic at the points A x = 0.5 
and A x — 1.5. These two points satisfy the condition 



J z 'Jx 



Jy and J x 



Jy + Jzi respectively. It is obvi- 



ous that these two points are on the critical lines denoted 
by the thick dotted and solid lines in Fig. [51 



III. ENTANGLEMENT 

It has been shown that the entanglement also exhibits 
critical behavior at the quantum phase transition point 
for both spin (see, e.g., [22L l23l . [2J|) and fermionic sys- 
tems (see, e.g., [H, HU). Also, it was shown [13, S| 
that there is a general relation between the bipartite en- 
tanglement and the quantum phase transition. In this 
section, we show that the nonanalyticity of the ground- 
state energy in the extended Kitaev model results from 
the correlation functions [see Eqs. ([24]), ([28]) and ([29]) 
for their definitions]. Furthermore, we show that the bi- 
partite entanglement also exhibits nonanalyticity at the 
quantum phase transition point and its nonanalyticity 
is also due to the nonanalyticity of the same correlation 
functions. This reveals that both the ground-state en- 
ergy and the bipartite entanglement can characterize the 
quantum phase transition in the Kitaev model. 
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1 + denotes (A^ — A y ) — > 1 with (A^ — A y ) > 1, and 1" 



A. Correlation functions and nonanalyticity of 
ground-state energy 

From Hellmann-Feynman theorem [29j , we have 
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where E is the ground-state energy per site given in 
Eq. (15), -ff is the Hamiltonian (5) (rescaled by J z ), Tr 
denotes the trace over the ground-state subspace, and 
P = 1.9) (ff I is the density matrix of the system. When 
|g)(<7| is traced over all spins except the two spins at r n . m 
and T n ',m't t ne reduced density matrix is 

p{r n , m) Fn',m') = T^'(\g)(g\) 
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where c" (er Q ) are Pauli matrices tr 21 , cr 2 ' and a z for a 
(a')= 1 to 3, and the unit matrix for a {a')— 0. When 
the two spins at r n m and r n+l rra are linked by an x-type 
bond, the reduced density matrix becomes 

/ \ ^ / I 21 X I \ X X 

P{r n ,rn, Yn+l,m) = — \g \o~ n ^ m o~ n+lm \ g) o~ nm o~ n+lm 
1 
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where n + m is an odd integer, and I is the unit operator. 

Because of translational invariance, the correlation 
function 

Qx = (S|<m<+l,m|ff) (24) 

is spatially invariant. Thus, Eq. (|23p can be written as 
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where n + m is an even integer for both Q y and £/ z . Here 
Eqs. (|25|) -([27 | are the results obtained for the reduced 
density matrix when the two spins at r„ iTO and r n / )m < 
are nearest- neighbors. When the two spins at r„ im and 
r n',m' are not nearest-neighbors, the density matrix is 
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Using the Jordan- Wigner transformation (3) and the 
definitions (D and (5) for the Majorana fermions, we can 
derive that 
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From Eqs. (TO), (12), and (3TJ we have 
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with -1 < Q y < 1. 

From Eqs. (2T), (23), and (28), it follows that 
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Equation (35) shows that the directional derivatives of 
the ground-state energy per site are determined by the 
correlation functions Q a (a — x, y) and their derivatives. 
Section II shows that dE/dl and d 2 E/dl 2 are continu- 
ous, while d 3 E/dl 3 is discontinuous on the critical line, 
e.g., A x = 1 + A y (i.e., the thick solid line in Fig. [2). 
Equation (35) reveals that the nonanalyticity of E on 
the critical line is due to the nonanalyticity of Q a . As 
shown in Appendix B, 
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where T is given in Eq. (fH)|) . Equation (|3"T)) sho-v 
the spin-spin correlation function C/ Q can signal th( 
turn phase transition, similar to the bond-bond c 
tion function in the original Kitaev model 21 1. 
Eqs. (J35J) and l[3T|). we have 
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which is the same as in Eq. (|18|) . This further 
that the nonanalyticity of the ground-state energy 
from the nonanalyticity of the correlation functio: 



B. Nonanalyticity of entanglement 

We now focus on the bipartite entanglement of the 
ground state \g) between two spins (at r„ jm and r n / >m ') 
and the rest of the spins in the system. We use the von 
Neumann entropy to measure the entanglement between 
these two spins and the rest of the spins in the system. 
The von Neumann entropy can be defined by [3fj| 

S a = -Tr [p(r n , m , r n+ i, m ) log 2 p(r n , m , r„+i >m )] , (39) 

where Tr denotes the trace over the two-spin Hilbert 
space, and a — x if n + m is an odd integer, and a = y 
if n + m is an even integer. Also, this entropy can be 
written as 



S a = - Xi log 2 Xi, 



(40) 



where the sum runs over the four eigenvalues Xi of the 
matrix p(r n , m , r n _|_i )m ). From Eqs. ([2^)1 and (|2"5|) , it fol- 
lows that 



Ai = a 2 = j(i + g a ), 
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FIG. 4: The bipartite entanglement (i.e., the von Neumann 
entropy) S x and its first, and second derivatives with respect 
to A x , where A y — 0.5, A = 0.1 and tp = (which corresponds 
to the horizontal thin dotted line in Fig. [2}. Obviously, S x 
and 8S x /dA x are continuous functions, but d 2 S x /dA x is dis- 
continuous at the transition points = 0.5 and 1.5. 



Thus, we have the entanglement measure 
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which is determined by the correlation function g a , sim- 
ilar to the thermal entanglement [3l| . 

To see the relationship between the entanglement and 
the quantum phase transition, we analyze the directional 
derivatives of the von Neumann entropy with respect to 
the driving parameters along any direction I. The first 
and second directional derivatives of the bipartite entan- 
glement are 
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From Eqs. [[36] ) . ([37 ]) . (|42 |) . and we have 
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and 



that this model displays a continuous quantum phase 
transition on the critical lines separating the Abclian and 
non-Abelian phases, where the third derivative of the 
ground-state energy is discontinuous. Also, we use the 
von Neumann entropy as a measure of bipartite entan- 
glement to study this topological quantum phase transi- 
tion. Our results show that the bipartite entanglement is 
also nonanalytic on the same critical lines as the ground- 
state energy. Moreover, we show that the discontinuity 
of the second derivative of the bipartite entanglement is 
related to the discontinuity of the third derivative of the 
ground-state energy. Our approach directly reveals that 
both the entanglement and the ground-state energy can 
be used to characterize the topological quantum phase 
transition in this Kitaev model. 



d 2 S x 


d 2 S x 




dP 


1+ dP 


1- 


d 2 Sy 


d 2 Sy 




dP 


1+ dp 


1- 




where 



Got ~ Ga | 2+ — ^ a 1 1" 



£,(46) 



(47) 



Equation (|4l)|) shows that the bipartite entanglement 
is nonanalytic with its second directional derivative 
d 2 S a /dP discontinuous at the critical line A x = 1 + A y 
(denoted by the thick solid line in Fig. Because 
— 1 < Qa < 1) it follows from Eq. (|43[) that the disconti- 
nuity of d 2 S a /dP is due to the discontinuity of d 2 Q a /dP. 
Similarly, it can be shown that d 2 S a /dP also exhibits 
a discontinuity on the critical lines A y = 1 + A^ and 
A x = 1 — A y (denoted by the thick dashed and dot- 
ted lines in Fig. ^ which is due to the discontinuity 
of d 2 G a /dP on these lines. As in Fig. [3l we choose 
A = 0.1, if = and A y = 0.5 as a typical example to 
show S x and its first and second derivatives with respect 
to A. x (see Fig. HJ) - It is clear that S x and its first deriva- 
tive dS x /dA x are continuous as a function of A^, but 
its second derivative d 2 S x /dA 2 is discontinuous at the 
quantum phase transition points A^ = 0.5 and A x = 1.5. 

As shown above, both the nonanalyticity of the 
ground-state energy and that of the bipartite entangle- 
ment are due to the nonanalyticity of the spin-spin cor- 
relation functions. This reveals that the ground-state en- 
ergy and the bipartite entanglement are closely related 
with each other, and both of them can be used to char- 
acterize the topological quantum phase transition in the 
extended Kitaev spin model. 
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APPENDIX A 

This appendix focuses on the analyticity of the first, 
second and third directional derivatives of the ground- 
state energy on the critical line denoted by the thick solid 
line in Fig. [5] From Eq. ([15)) , the derivative of the ground- 
state energy with respect to A x is 



dE 
dA^ 



V3_ 
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■ £ (k) 



d 2 k- 



A 



BZ 



d2k a 
d e(k) 



e(k) 

A y/3 



16tt 2 



BZ-D ^(k) 



(Al) 



IV. CONCLUSION 

In conclusion, we have studied the topological quan- 
tum phase transition between Abelian and non-Abelian 
phases in the extended Kitaev spin model on a honey- 
comb lattice. From the ground-state energy, we show 



where 



A = A : , 



A y cos k x — cos 



k x ^/3ky 



(A2) 



and D denotes two small regions in the first Brillioun 
zone, i.e., half of the disk with radius e, which is centered 



8 



at (tt, — 7r/v3)) and half of the disk with radius e, which 
is centered at (—tt, 7r/\/3), where e«l. When A x — A y 
1 . e(k) becomes zero only at the points (=p tt, ±7t/\/3), so 
Jbz-d ^ 2 ^e7ET * s analytic because BZ — D is the region 
excluding D in the first Brillioun zone. For the integral 
in the region D, we can approximate it as 



d2k ^A 
d e(k) 



where 
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Similarly, 



as 
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(A8) 
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From Eqs. (16]), (|A7)) , and |[X8)) . it follows that 

_ dE 

i+ 



dE 
~dl 



dl 



(A9) 



Using the same procedure as above, we can obtain 
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where 
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Q1+Q2 sin(20 + 



(All) 



A = A* -Ay -I, 

51 = P 1 + P 2 sm(26 + <t> 2 ), 

5 2 = Qi + Q2sin(26» + 0i), 



(A4) 



and 

Qi 



A X Ay 



4A 2 + A, 



A„ 



1 



Q2 = \^{2A X A V + 8A 2 - A x + A. y ) 2 + i{A x + A y f 
Pi 

P2 = WA 2 -A, + f. 

' 2A x A y + 8A 2 - A x + A y 
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!>i = arctan 
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From Eq. 



V3 
3l, we have 
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where 1 + denotes A x — A 
A„ 



1 with A T , — A y > 1, and 



1 with A 7 



A. y < 1. When 



1 denotes A^ 
e -> 0, it follows from Eq. (JXSJ that / BZ _ D d 2 ^^ -> 
on the critical line A x = 1 + Ay (i.e., the thick solid line 
in Fig. [J). Thus, from Eq. (|A1[) . we have 
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APPENDIX B 

This appendix gives results regarding the analyticity of 
the correlation function Q a (a — x, y) and its first, and 
second directional derivatives on the critical line denoted 
by the thick solid line in Fig. fj] Similar to Eq. (|A1|) , one 
can divide the integral in (|32[) into two parts: 

Qx = -^-o / a k—— 
8tt 2 Jbz e(k) 



V3 
~8tt 2 



d 2 fc 



A V3 
e(k) 8vr 2 



BZ-D e l K J 



,(B1) 



where A is given in Eq. (|A2[) . Using the same procedure 
for Eqs. (|A9p and (|A10[) . we can derive from Eq. (|B1|) 
that 
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where T is given in Eq. ()A11[) . 
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From Eq. (|34|) . we can obtain 
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